Exact dynamics for fully connected nonlinear networks 
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We investigate the dynamics of the discrete nonlinear Schrodinger equation in fully connected networks. For 
a localized initial condition the exact solution shows the existence of two dynamical transitions as a function of 
the nonlinearity parameter, a hyperbolic and a trigonometric one. In the latter the network behaves exactly as 
the corresponding linear one but with a renormalized frequency. 
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I. INTRODUCTION 

Nonlinear dynamics in complex networks incorporates 
competition of propagation, nonlinearity and bond disorder 
and may find some applications in complex natural or man- 
made materials. For a given random network of N sites the 
limit of fully coupled lattice where each site is connected to 
every other one with the same strength plays an important role 
since, for zero nonlinearity, introduces a high degree of de- 
generacy and thus localization |l|]. It is therefore interesting 
to probe the dynamics in this fully coupled limit when the net- 
work is nonlinear. We use the Discrete Nonlinear Schrodinger 
(DNLS) equation for this study since the latter is a prototypi- 
cal equation with a large number of applications in condensed 
matter physics, optics, Bose-Einstein condensation, etc. |0-Bl- 

The DNLS equation in the fully coupled or Mean Field 
(MF) limit is given by: 
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where \|/„(f ) is a complex amplitude, V denotes the constant 
overlap integral connecting any two sites, e„ a local site value 
and y is a parameter that controls the strength of the nonlinear 
term. When t is time, the equation describes the dynamics of 
interacting particles with the nonlinear term representing po- 
laronic effects or interaction of bosons in an optical lattice. 
When t is a space variable, the equation describes electric 
wave modes propagating in fibers that include nonlinearity. 
In the present analysis we take e„ = while the DNLS norm 
equal to one, viz. L|V«| 2 = 1 ■ F° r notational simplicity we 
rescale time to x = Vt and y to % = y/V; we note that for % > 
we have the defocusing DNLS equation while for % < we 
obtain the focusing case. 



out that no selftrapping occurs for positive nonlinearity val- 
ues. This dynamical behaviour is analogous, but not identical, 
to the well known selftrapping transition of the dimer JH . 
Andersen and Kenkre [8|] solved analytically the symmetric 
trimer case for localized as well as symmetrically delocalized 
initial conditions. They wrote explicitly the time evolution 
at selftrating (% = —6), connected it to the stationary states, 
calculated the averaged occupation probability of the initially 
populated site and showed that for a given value of the non- 
linearity parameter extreme derealization occurs. They also 
pointed out that the nature of selftrapping in the trimer was 
distinct from that in the dimer case. Subsequently, Andersen 
and Kenkre analysed the fully connected A^-site problem 
using similar methods for localized and partially localized ini- 
tial conditions and showed that a "migratory" transition takes 
place in some cases. The latter is related to the fact that local- 
ization of the initially non-occupied sites may occur for some 
parameter values. In |§ a general expression in terms of 
the Weierstrass elliptic function is given for the time evolu- 
tion of the occupation probability as well as expressions for 
some special cases. Additionally the time-averaged occupa- 
tion probability is analysed. Finally, Molina ITIOII used an 
extension of the method of ref. Q] and determined explicitly 
the values of the nonlinearity parameter at which there is a 
"transition" for arbitrary values of N. 

From the previous review of the existing literature it is clear 
that the problem of nonlinearity induced localization in the 
fully coupled, MF limit, has been analysed thoroughly and 
many of the features of the solution have been already discov- 
ered. However, we believe that a complete dynamical picture 
has not emerged yet. Thus we will follow the approach of 
ref. |@], re-derive the complete dynamical solution and use it 
to focus on the explicit dynamics and not just to the average 
occupation probability. 



II. REVIEW OF EARLIER LITERATURE 

The problem of the DNLS equation in fully coupled lat- 
tices has been addressed in the past JH p7l- Hoil . Eilbeck et 
al. analysed the stationary states of fheproblem for the 
general A^-site case. Molina and Tsironis J7[] showed that in 
the case of N = 3, i.e. a symmetric trimer, and for a local- 
ized initial condition, a transition occurs for % = —6 where 
the evolution ceases to be periodic and the averaged proba- 
bility becomes equal in all three sites. It was also pointed 



III. NONLINEAR LOCALIZATION 

In the purely linear case for % = 0, Eq. (Q~|) reduces to a lin- 
ear eigenvalue problem that has N — 1 fold degenerate eigen- 
values This high degree of degeneracy forces an initially 
spatially localized state to remain localized executing incom- 
plete oscillations to the other, initially unoccupied, sites. This 
linear localization at the initially excited site increases with 
the system size. 
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A. Nonlinear mean field limit 

For the localized initial condition we observe that if we set 
initially all probability at site one, then due to the symmetry 
of the lattice, the subsequent evolution of all sites except the 
initially populated one will be identical; as a result we can cast 
the original Adequation problem to a simpler one consisting of 
only two equations, viz. I7 H10I1 

iVfi = (N- l)y 2 -xlVi|Vi (2) 
iVf2 = ¥1 +(A ? -2)i|/ 2 -xl¥2| 2 ¥2 (3) 



with \\\t\ \ 2 + (N — 1)|¥2| 2 = 1- To obtain the set of equations 
[3] we have set \|/2 = ¥3 — •■ = ¥w at a U times. We proceed 
by introducing the variable p(t) = |\|/i| 2 — (N — 1)|¥2| 2 and 
cast the set of complex equations (Q]) first into a second order 
equation that has the form of Newton's law and subsequently 
into an equivalent potential problem for the localized initial 
condition p(0) = 1 JH3-0, viz. 



- U(p) = a p 4 +4aip 3 + 6a 2 p 2 + 4a 3 p + a 4 (4) 
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We now form the invariants g2, g3'. 

g2 = ao«4 — 4flifl3 + 3a 2 (12) 

g3 = a§a 2 a\ + 2aifl2fl3 — a\ — aoa 2 — a\a 4 (13) 

The solution is given in terms of the Weierstrass elliptic func- 
tion T (t , g2 , gi ) as follows UlTl : 

24(N- 1) 

P(T) ~ 1 ~ 12fP(T,g 2 ^3)+X 2 + 2(^V-2)x+^ (14) 



This solution depends only on the nonlinearity parameter % 
and the number of sites N and presents the exact solution of 
the nonlinear DNLS problem in the MF limit for localized 
initial condition and arbitrary site number; it is equivalent to 
the solution in ref. {gj for perfect initial localization on a 
single site. The Weierstrass function is doubly periodic with 
periods 2co and 2co' where co/co' is not real. The discriminant 
A = g\ — 21 g\ may be expressed as follows 



A(N) = ~^-N 2 X 2 g(N)( X + 2N) (15) 
to 

g(N) =2(N- 2)x 3 + 2 (N - 2) (3N - 8) % 2 + (6N 3 - 53N 2 + 1 12N - 64) % + 2{N - 1 )N 3 (16) 
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When the discriminant becomes zero, the Weierstrass solu- denoted hereafter as %o, % c and %[. for %q = 0, g(% c ) = and 
tion changes character signalling a change in the physical be- yj c = —2N respectively, 
haviour of the system. The quantity A(N) has three real roots, 
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The root for %o = corresponds to the linear MF prob- 
lem resulting in g2 = N 4 /12, gj, = N 6 /216, the Weierstrass 
function becomes a trigonometric function, viz. £P (t,g2,g3) = 
-N 2 /12 + N 2 /4sin 2 [(Nt/2)] and the solution for p(x) is 



p(%) =p,in{x) 



{N-2) 2 +4(N-\)cos(Nx) 

n 2 



(17) 



a result that agrees with the linear solutions derived in ref . JTj] . 
We note that in the linear case the system executes oscillations 
between p = 1 and p, inM = [{N -2) 2 -4(N - 1)] /N 2 while 
p(%) oscillates around < pu„ >= (1 - 2/N) 2 . Only for N = 2 
we have complete oscillations in the p-range 1 , — 1 while for 
all other geometries the linear system oscillations are incom- 
plete, viz. there is always some occupation probability at the 
originally populated site. In the limit of very large N there is 
extreme localization at the initially populated site. 

In order to understand the role of % c , %' c in the dynamical 
evolution, it is instructive to address first explicitly the regime 
for N < 4; although these cases have been studied in the past, 
the present exposition furnishes additional results and new in- 
tuition. 



B. 2<N<4 

In this range we have %' c < % c with the equality holding for 
N = 4. 



1. N=2 

In the dimer case we have % c = 4 and %' c = —4; when % 
assumes one of these two critical values, the invariants be- 
come g2 = 4/3, gi = —8/27 , the frequencies CO = K(l) = °°, 
co' = ;K'(1) = in/2 and the Weierstrass function becomes a 
hyperbolic function, viz. 2>(f, 4/3, -8/27) = 1/3 + l/sinh 2 f 
leading to p(t) = sech(2t). This situation corresponds to the 
well known selftrapping transition for a localized initial con- 
dition where the solution at the transition is expressed as hy- 
perbolic secant @]. We note that in the dimer case % c = \%' c \ 
and for — % c < %< % c the particle executes complete oscilla- 
tions. 



2. N=3 

In ref. |01 a transition was found occurring at %' c = —6, 
while in refs. |0| HI] the dynamical study was performed. 
Analysis of the discriminant shows that in addition to the pre- 
viously identified root %' c = —2N = —6 there is also the real 
rootofg(x c ) =0at 
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For % = Xc = — 6 the invariants are g2 — 3/4, gi 
1 8. co becomes infinite and the Weierstrass function turns 
into a hyperbolic function, viz. 3>(/,3/4, — 1/8) = 1/4 + 



3/4 sinh( y 



The time evolution for p(x) be- 



3-4im/7 2 ( v /374T) 
1 +4sinh 2 ( x /J/4x) 



(19) 



At long times p{x)u mx ^ x = — 1/3 rj —0.333 corresponding 
to equal probability 1/3 on all three sites and should be com- 
pared with pim,M = —7/9 ps —0.777 and < pu„ >= 1/9 ps 
.111 When x = %c both invariants are positive and the Weier- 
strass function turns into a simple trigonometric function. 

We now describe the dynamics of the trimer as a function 
of the nonlinearity parameter %. For x = 0, p{%) performs in- 
complete oscillations between p(0) = 1 and ph„.m = —7/9. 
Increasing % to positive values augments gradually the effect 
of localization and the particle stays almost completely at the 
initial site for large values of the nonlinearity parameter. For 
negative values of %, however, the behaviour is markedly dif- 
ferent; upon increasing % in absolute value, the particle be- 
comes more detrapped and performs oscillations beyond the 
piin,M limit. This behaviour continues as we further increase 
|X| and at x = — 4 the particle executes complete oscillations 
to the other sites. Further increase of nonlinearity reduces the 
amplitude of the oscillatory motion while the oscillation pe- 
riod increases. For % = %' c there is a transition to a hyperbolic, 
non oscillatory behaviour; the particle collapses from com- 
plete oscillations to complete equipartition of probability in 
the three sites (at infinite times). This is a remarkable transi- 
tion since it localizes again the particle and all motion stops 
asymptotically 0Htl- As the value of |x| increases further, the 
trapping tendency increases rapidly and the particle becomes 
more trapped in the initial site. 

As the trimer crosses the value % c the discriminant becomes 
zero again with positive invariants; this means that the time 
evolution simplifies to one with trigonometric evolution. Al- 
though this feature is not exceptional in the trimer, it is quite 
interesting when it appears for > 4; we will comment on 
this trigonometric transition below. 

The transition at %' c is a selftrapping transition in the sense 
that it forces probability equipartition on all three sites, similar 
to the corresponding behaviour in the dimer. It is remarkable 
that nonlinearity first detraps the particle and after reaching a 
regime of complete oscillatory motion it then forces the par- 
ticle to selftrapping. As a result nonlinearity first restores the 
degeneracy that has been lifted by the presence of more than 
two sites and subsequently acts in way analogous to the dimer 
case, viz. reaches an equipartition state through a hyperbolic 
function transition. 



3. N=4 



The tetrahedral configuration obtainned for = 4 is rather 
interesting with % c = %' c = —8, g2 = g3 = and the Weierstrass 



4 



function becoming a simple rational function T (t , 0, 0) = 1 /t 2 
leading to the solution 



/»(*) 



l-2x 2 
1+4t 2 



(20) 



For long times the solution reaches the value —1/2 while 
Piin.M = — 1/2 and < /;>/,■„ >= 1/4 . For positive nonlinear- 
ities the localization tendency increases while for negative % 
we have first detrapping, complete oscillations for % = —6 and 
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subsequent asympotic equipartition at % = —8 with the transi- 
tion being now algebraic as seen in Eq. d20"l i. Further absolute 
increase of % produces nonlinear localization. We note that 
as in the trimer also in the present case at % = %' c we have an 
asymptotic equipartition of the probability on all four sites; in 
that sence we could use the term "selftrapping" meaning how- 
ever the asymptotic and irreversible onset of equal probability 
on all sites. 
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C. yv>4 

We observed from the previous three cases that in the dimer 
the two roots at % c and %' c determined two distinct, yet iden- 
tical in nature, selftrapping transitions. In the trimer yj c was 
responsible for the hyperbolic transition, % c for trigonomet- 
ric evolution, while the fact that the two roots were equal in 
the tetramer turned the (unique) transition to an algebraic one. 
The (negative) roots % c and %' c "collide" in the tetramer and 
we observe that their role changes for N > 4. Specifically, we 
find that \%' c \ > \% c \, with the algebraic transition occurring at 
% c while at %' c we have a specific trigonometric evolution. 

For negative nonlinearities as we increase in absolute value 
X we have initially detrapping, complete oscillatory motion 
followed by subsequent reduction of the amplitude of the mo- 
tion that leads to the hyperbolic transition at % c . This state 
does not populate equally at long times the sites of the A^-mere 
as in the cases for N = 2, 3, 4; the asymptotic population of 
the initially excited site is now larger than that of the other 
sites. Further increase of |%| leads to periodic motion with in- 
creasing amplitude reduction. As the nonlinearity parameter 
crosses the value %' c = ~2N the dynamics becomes trigono- 
metric while, at the same time the amplitude of oscillation 
becomes identical to that of the corresponding linear A^-mere. 
Thus at %' c the nonlinear system behaves identically to the cor- 
responding linear one in what regards the site probabilities, 
nevertheless the evolution is done with different frequency. 
Further increase of nonlinearity leads to more localization. 

The value of % c is determined numerically from the solu- 
tion of the cubic equation g(% c ) = 0; the values obatained 
coinside with those found in ref. IToll . We observe that as N 
increases this root separates from the one at —2N and asymp- 
totically reaches the value — N. Thus in a very large system 
we expect to have the hyperbolic selftraping at nonlinearity 



values of order —N while at —2N the nonlinear system be- 
haves as the corresponding linear one with renormalized fre- 
quency. We also observe that the system reaches a regime of 
complete oscillations for %' = —2(N— 1) while %' < %' c <% c \. 
For N > 15 the equality for %' does not hold any more and 
the system begins to execute incomplete oscillations reaching 
a maximum at %' < % £ . Thus, while for N < 14 the parti- 
cle first reaches complete derealization and at a later stage 
relaxes hyperbolicaly, for N > 15 the two phenomena occur 
almost for the same %-value, first reaching a maximum in the 
derealization and subsequently, for infinitesimal change in 
the x-value, hyperbolic -type evolution. In this last regime the 
maximum excursing occurs for % very close to the hyperbolic 
transition value % c . 

The exact expression for the evolution at the trigonometric 
transition at %' c = —2N for N > 4 is given by: 



4-6N + N 2 + 2{N-2)cos(^/N(N-4)x 



N 



N-2-2cos[^N(N-4)z 



(21) 



The oscillation extends from p max = 1 to p min = (A^ 2 — 8A^+ 
S)/N 2 ; these values are identical to the one obtainned for the 
linear MF limit. Thus, in the MF limit of the DNLS equa- 
tion there are two values of the nonlinearity parameter, viz. 
X = and at which an initially localized excitation os- 
cillates trigonometric ally with the same amplitude, but with 
completely differenent frequencies. In the second case, the fi- 
nite value of nonlinearity is equivalent to a "polaronic" dress- 
ing of the excitation leading to a smaller oscillation frequency 
compared to the purely linear oscillation at % = 0. 

The hyperbolic transition at % c is given, on the other hand, 

by 



p(x) = 1 



24(AT-1) 



12c 



sinh 2 [-v/3cx] 



-X C 2 +2(N-2)% C +N 2 



(22) 



where c is the value of the double root of the discriminant of the equation x 3 — g2X — gi — while the other root is equal 
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more difficult. The time evolution for various networks is pre- 
sented in Fig. (1) while the special nonlinearity values % c and 
Xc as a function of the system size N are portrayed in Fig. (2). 



FIG. 1: (Color online) Exact dynamics for fully connected networks: 
probability p(%) of Eq. d 1 4b as a function of the rescaled time x. Each 
subfigure is for a network with different size and specific nonlinear- 
ity. In all plots the periodic continuous red line presents the motion 
for localized initial condition and for the specific N and % while the 
dashed orange line the solution of the corresponding linear network. 
The continous green line is the solution at % = %' c = —2N while the 
dashed blue line the solution at % — % c . The three horizontal lines 
parallel to the time axis denote the maximum excursion for the cor- 
responding linear network pu„,m (blue dashed), the equipartition line 
at 1 — 1 /N ( purple) and the maximum range for the probability p(x) 
at —1 (dashed with spaces). In (a) N = 3 and x = — 4. We note 
the restauration of complete oscillations to —1 due to nonlinearity 
and the trigonometric evolution (blue curve) that in the present case 
is distinct from the purely linear motion (dashed orange curve), (b) 
N = 4 and x = — 6. The green and blue curves collapse onto a sin- 
gle algebraically decaying solution (% c = x' c = ~ 8 that relaxes to the 
equipartition line, while complete periodic motion to the other sites 
is also observed (red curve reaching — 1). In this case the maximum 
linear excursion also reaches to the equipartition line, (c) N = 10 and 
% = -18. The solutions for the roots at Xc - - 19.333 and %' c = -20 
have now changed roles with the former (blue line) designating the 
hyperbolic transition (to non-equipartition- equipartition is denoted 
by the dashed purple line) while the latter (green line) forcing the 
network to pure trigonometric evolution with the same amplitude as 
for the linear (% = 0) case. In (d) N = 15 and % = — 15 we observe 
similarly the hyperbolic trannsition ( Xc — —27.9772, blue line) to 
non-equipartition as well as the two trigonometic evolutions for x = 
(dashed orange line) and % r c = — 20 (green line) 



to —2c. 

The presence of focusing nonlinearity induces complete or 
partial detrapping to the MF network followed by retrapping 
for larger |%|. The exact solution of the nonlinear model for 
arbitrary N shows thus very interesting features of detrapping, 
hyperbolic relaxation as well as effective linearization of the 
nonlinear system. As the network size N grows, these interest- 
ing fearures appear for increasingly larger values of the non- 
linearity parameter |%|, a feature that makes their observation 
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FIG. 2: (Color online) Plot of the values — Xc (red) and — %' c (green) 
as a function of the network size for two regions: Left panel, for 
N < 12 and right panel for 50 < N < 60. The black circle designates 
the "collision" of the two branches at N = 4. For large-/V Xc ~ — N 
IToll . while x'c = — 2# for any system size. 



IV. CONCLUSIONS 

The analysis of the MF model for the DNLS equation with 
a localized initial condition and arbitrary number of sites 
shows that the system is rich and has interesting dynamical 
behaviour. In the linear case (% = 0) the presence of the long 
range bonds leads to localization that increases with Af JH]. 
As focusing nonlinearity is turned on we observe initially a 
tendency for detrapping, i.e. the particle executes oscillations 
of larger amplitude that can reach complete or partial (N > 1 5) 
escape from the initially populated site. 

The time dependence of the system is typical of elliptic 
function dependence, viz. periodic oscillatory motion that 
includes an infinite number of frequencies. This time de- 
pendence changes for two specific values of the nonlinear- 
ity parameter at % c and %' c respectively. For the former value 
the time evolution becomes hyperbolic (except in the trimer 
where the hyperbolic evolution occurs at %' c ) and periodic mo- 
tion ceases at long times. In the well studied case of N = 2, 
the state for this critical nolinearity separates periodic delo- 
calized from localized motion and we have selftrapping. For 
N = 3, 4, while there is no such separation in the hyperbolic 
transition, there is asymptotic equipartition of probability for 
this value. For N > 4 the hyperbolic solution at % c simply rep- 
resents an asymptotic state where motion ceases while the ini- 
tially populated site retains most of the probability, a tendency 
that increases with N. Thus, if by "selftrapping" we mean 
the separation of delocalized to localized motion then, strictly 
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speaking, it only occurs in the dimer case. On the other hand 
we may extend the term to include the case where asymptotic 
equipartition of probability takes place; in this case the sym- 
metric trimer and tetramer also qualify, but not networks with 
larger number of sites. 

A new feature appears at % = %' c (% = % c in the trimer); the 
elliptic function evolution becomes trigonometric and the in- 
finite set of frequencies it encompasses collapses to a single 
one. Additionally, for N > A, the amplitude of the oscillation 
becomes identical to the corresponding one for % = 0. Thus 
the evolution at %' c is effectively linear but with a new fre- 
quency equal to y/N(N — 4) (for N > 4) while in the purely 
linear case the frequency is N d]. If we think that the nonlin- 
ear term stems from the coupling with additional degrees of 
freedom lfl2ll then the evolution becomes more sluggish since 
the particle is dressed and carries with it the additional degrees 
of freedom. It is remarkable that such a dressing leading to 



frequency renormalization may occur and still the evolution 
be trigonometric with the same amplitude as in the purely lin- 
ear case. In the large-Af limit the two frequencies approach 
one the other. 

As the size of the system increases we observe two distinct 
features, one at % c ~ —A? where time-hyperbolic behaviour 
takes place while at %' c = —2N where we have an effective 
linearization of the motion. Both features appear after the ini- 
tial detraping and retraping has occurred as a function of |%|. 
The onset of nonlinear localization is typically associated with 
discrete breather solutions that are localized in space and time 
periodic. The case for % = yj c found here would then corre- 
spond to a "linear breather", i.e. a localized nonlinear solution 
with exact trigonometric evolution. It would be interesting to 
find out if a similar situation can occur in other discrete non- 
linear systems as well. 
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